Inter-layer Hall effect in double quantum wells subject to in-plane magnetic fields 
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We report on a theoretical study of the transport properties of two coupled two-dimensional electron 
systems subject to in-plane magnetic fields. The charge redistribution in double wells induced by 
the Lorenz force in crossed electric and magnetic fields has been studied. We have found that 
the redistribution of the charge and the related inter-layer Hall effect originate in the chirality of 
diamagnetic currents and give a substantial contribution to the conductivity. 
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I. INTRODUCTION 

Semiconductor quantum wells, multi-wells and super- 
lattices are designed to have an artificial electronic struc- 
ture determined by their construction rather thanJpy the 
properties of individual semiconductor materials p. The 
band profiles of multi-wells are formed by a sequence of 
quantum wells separated by barriers. Electrons condense 
to two-dimensional (2D) electron layers localized in the 
wells and the tunnel-coupling of the layers leads to the 
formation of new electron systems, rich in optical and 
transport effects, that can be modified by application of 
the magnetic field. 

Two experimental arrangements of the magnetotrans- 
port measurement in multi-well structures subjected to 
the in-plane magnetic field B are possible: with the cur- 
rent /flowing perpendicular to the magnetic field (J* _L B) 
or parallel with it (J* || B). The Hall effect plays a role 
in the case J_L B. fn wide macroscopic samples (super- 
lattices) the standard method of the description of the 
Hall effect, via the inversion of the conductivity tensor, 
applies. This is no longer true for the double- well struc- 
tures the width of which is microscopic. This is why we 
present here a novel theoretical approach to the calcu- 
lation of the Hall effect, based on its textbook defini- 
tion: the Hall electric field is caused by the accumulation 
of positive and negative charges on the opposite sample 
edges, induced by the Lorentz force. 

The difference between conductivities along and per- 
pendicular to the magnetic field direction was investi- 
gated experimentally in several papers, □□. rThe results, 
obtained by solving the Boltzmann equation have been 
found in a qualitative agreement with the experimental 
observation. However, no attention has been paid to the 
effect of the sample polarization due to the transfer of 
electrons between wells - to the Hall effect. 

Let us consider the setup sketched in Fig. 13. The elec- 
trons are pushed by a Lorenz force from the left to the 
right well. Resulting non-equilibrium charge distribution 
gives rise to a potential that compensates for the Lorenz 
force in a similar way as the Hall voltage does in macro- 
scopic samples. In a limiting case of very narrow wells, 



the double-well system can be viewed as a parallel plate 
capacitor: the induced potential reduces to an inter- well 
Hall voltage Uu, a potential difference between electron 
systems in individual wells. 




FIG. 1. Schematic picture of a double-layer system. Direc- 
tions of applied electric and magnetic fields are indicated. 

Note that in such an experimental arrangement the 
electron transport is closely related to the chirality of 
diamagnetic currents flowing in opposite directions in the 
left and right wells. In the thermodynamic equilibrium, 
the two currents are exactly balanced giving a zero net 
flow of electrons through the sample. In non-equilibrium 
states the electron system is polarized, more electrons are 
in the right well and an incomplete compensation of dia- 
magnetic currents results in a finite transport current. 
The aim of this paper is to show that this usually ne- 
glected correction to the in-plane electrical conductivity 
gives a significant contribution for realistic double well 
structures. 

For the sake of simplicity werfimploy the model of a 
double well described in a paper Two coupled, strictly 
two-dimensional electron layers confined in very narrow 
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potential wells are considered. Short-range scatterers, 
randomly distributed in individual wells, are assumed to 
be responsible for a finite relaxation time of electrons. 

We use the self-consistent linear response theory (Kubo 
formula) in the random phase approximation (RPA) (see 
e.g. u) to obtain the non-equilibrium charge distribution, 
the inter-layer Hall voltage and the conductivity tensor 
components. 



II. MODEL HAMILTONIAN 

In our model, the vector potential A = (zB, 0,0) is 
used to describe the influence of the in-plane magnetic 
field B — (0, B, 0) on the electron structure of a double- 
layer system. The corresponding Hamiltonian reads 



Hn = 



1 



1 



2m* 



(p x + \e\B Z y + —(f y +p 2 z )+V(z), (1) 



where m* denotes the electron effective mass and V(z) is 
a confining potential with two deep minima at z = id/2, 
V(d/2) = V{-d/2). The subscript denotes that no 
scattering mechanism is included. The form of (Q) al- 
lows separation of variables, the x, y dependent part of 
eigenstates of Hq is the plain wave with the two- 
dimensional wave vector k — (k x , k y )]^ To describe elec- 
tron eigenstates depending on the remaining third coor- 
dinate z we restrict ourselves to the tight-binding approx- 
imation employirui|-pnly the lowest eigenstates \tp a (z)} of 
uncoupled wells Uu. The index a — R, L distinguishes 
the right and left well. Thus, a complete single-layer 
eigenstate reads \ip a , k) = \ip a (z)) \k) and 



Pa = ^ \ fa,k)(ip a ,k\ 



(2) 



has a meaning of a projection operator to the states in 
the well a, Pl+Pr = 1- The Hamiltonian ([!]) is diagonal 
in the index k in our restricted basis and takes the matrix 
form 



where E^^^k) are single-well eigenenergies, 

2 *2, 2 



E L , R (k) 



If I V h k± 



(3) 



(4) 



The hopping integral t is given by t = {tp l\V (z)\tp r) = 
(ipit\V(z)\(pL)- Note that we neglected the effect of the 



*In the rest of the paper all the vectorial notation corre- 
sponds to two-dimensional vectors lying in the layer planes. 
The only exception are the vectors A and B introduced above. 



in-plane field on the energy spectra of individual layers 
and that the z coordinate has been replaced in (H) by 
±d/2 for the right and left layer, respectively. 

The diagonalization of the matrix (J3j> yields a pair of 

new eigenstates \x^t\ z ))- They depend only on the x 
component of k and the index i equals b for the bonding 
state and a for the antibonding state. The new basis 
\Xi^) is related to the basis of single-well eigenstates \tp a ) 
by 



-^VT±a\<pl)±^VTtA\<pr), (5) 



where A = S/y/8 2 + 1 2 and 5 is the abbreviation for 
h\e\dBk x /2m* . Then the eigenenergies of (|l|) reads 

h 2 kl 

(6) 



E b , a (k) = E$(k x ) + 

h 2 kl 



2m* ' 
e 2 d 2 B 2 



T VS 2 +t 2 . 



2m* 8m* 

The model is characterized by two parameters, the 
inter-layer distance d and the hopping integral t. It is a 
good approximation when the cyclotron energy and the 
tunnel coupling are small in comparison with the single- 
well quantization energies. 

The Hamiltonian H of the system with short-range 
scatterers is obtained by adding the impurity potential 
Vi m p to Hq. As mentioned above, we assume that impu- 
rities are distributed in both the left and right wells at 
random. To be more specific, the scattering on an indi- 
vidual impurity is considered to be intra-well (diagonal 
in the layer index a) and isotropic in k direction. The 
finite lifetime and transport relaxation time result from 
replacing the quantities characteristic for a given config- 
uration of scatterers by these quantities averaged over all 
possible configurations. 



III. CHARGE REDISTRIBUTION 

Let us first turn attention to the self-consistent pro- 
cedure of establishing the non-equilibrium charge distri- 
bution and the Hall potential Uh in the RPA, as this is 
the novel feature of our approach to the theoretical de- 
scription of magnetotransport in microscopically narrow 
systems. The basic scheme of our derivation is outlined 
below, more details are presented in the Appendix. 

The standard Kubo formula describing the linear re- 
sponse to a time-dependent perturbation operator H p is 
employed, our formalism is close to that used intl In our 
case, the operator H p is given as a sum of two parts: the 

externally applied potential represented by Hp 1 ^ and the 
self-consistent electric field of non-equilibrium electrons 
described by Hp 2 \ 

The external homogeneous electric field £ = (£ x , £ y ) 
is included into the vector potential A. Then the corre- 
sponding part of Hp takes the form 
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— v ■ 



£e l 



(7) 



where v — (v x , v y ) is the velocity operator and uj denotes 
the frequency. This perturbation yields the static polar- 
ization of the sample and induces the stationary current 
in the limit u> — > 0. 

The second part of the perturbation Hamiltonian, de- 

noted as H p , describes the Hartree potential derived 
from the non-equilibrium electron density by solving the 
Poisson equation. In our model, this potential reduces 
to the potential difference Uh between two 2D electron 

~ (2) 

layers and the operator H p can be written as 



(P R -P L ) 



(8) 



Here z stands for an operator of the z coordinate in the 
representation of single-layer eigenstates \(p a ) and Uh has 
to be determined self-consistently. 

We introduce SQ a , the non-equilibrium charge density 
per unit area of a well a, as an expectation value of the 
operator 



= eP a = e 



\(p a , k){(f a , k\ 



(9) 



evaluated using the non-equilibrium density matrix. Due 
to the charge conservation 5Qr is equal to —SQl, and 
we can write SQ instead of SQr {—SQ instead of SQl), 
for convenience. 

In our model, the Hall potential calculated in the 
Hartree approximation, Uh, is related to the local non- 
equilibrium charge density SQ by a particularly simple 
formula 



SQ = e—. 



(10) 



This expression corresponds to the parallel plate capac- 
itor in which the potential difference between plates is 
fully determined by their distance, the excess charge on 
one particular plate and the dielectric constant e of the 
insulator between plates. 



Since H p is a sum of Hp 1 ^ and Hp"' , also the charge 
density SQ a can be written as a sum of two contributions 



(2) 



6Q& and SQ {2) . 



We get the response to Hp in the form 



SQW =ihe 2 £ x J dEf Q (E) 
x Tr( S(E - H) 



(11) 



pdG^El dG-{E) p 



dE 



dE 



in the limit u> 0. Here (• • •) denotes the configuration 
averaging and fo stands for the Fermi-Dirac distribution 
function. As expected, the applied electric field £ y does 
not induce any excess charge. 



The response to the operator H p 
tion (§fi, yields 



(2) 



defined by equa- 



SQ { <? =-^U H J dEf (E) 



Tr (S(E-H) zG-{E)P a +P a G + {E)i 



(12) 



The equations (Jl_lJ) and ([L2|) relate the induced charge 
densities to the electric fields £ x and —Unjd and their 
structure reminds the equation (|Io|). Therefore, if we 
consider them as definitions of the generalized dielectric 
functions e' 1 ' and e^ 2 \ we can write 



SQ 



(i) 



e (1) £ x , SQ 



(2) 



r(2) 



Uh 
d 



(13) 



where again the layer index is omitted and we have in 
mind the partial densities in the right well. 

Making use of SQ = 8Q^ + SQ^ 2 \ the decomposition 
of the local charge density SQ into two components, the 
self-consistent equation (|TJj) takes the form 



Uh_ 
d 



and its solution reads 

Uh_ = 

d £ (2) 



U H 



(14) 



(15) 



This formula relates the self-consistent Hartree field of 
electrons to the x component of the externally applied 
electric field £, as we have anticipated in the introduc- 
tion. 



IV. IN-PLANE CONDUCTIVITY 

With the established relations between Uh, SQ and 
£ x , the calculation of the conductivity is a straightfor- 
ward procedure. The components of the tensor a are de- 
termined through the expectation values of the current 
operator 



j = ev + 



e 2 £ 



(16) 



calculated using the non-equilibrium density matrix 
given by a response to the perturbation H p in the limit 
w->0. 

Similarly as the generalized dielectric function e, the 
conductivity tensor can be written as a sum of two contri- 

(i) 



n The response to Hp leads to the 



butions o^ 1 ) and V 2 ' 

standard expressions a for the static conductivity a^ 1 ). 
We can write for its components 



tQ> = irhe 2 dE - 



dfa(E) 
dE 



Tr (5{E - H)v s 5{E - H)v t 



(17) 



where s stands for x or y. 
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The response to the self-consistent field described by 

~ (2) 

the perturbation Hamiltonian Hp gives a non-zero value 
only for the current flowing along the x direction. The 
corresponding contribution to the conductivity reads 



r( 2 ) 



,(2) 



dEfo(E) 



(18) 



xTr (S(E-H) zG~(E)v x +v x G + (E)z 



Thus, while the conductivity a yy is determined only by 



the external field £ y and a yy 



a yy , in the case of the 



x direction also the self-consistent Hall field Uh helps to 
conduct the current and a xx 



u XX 



a {2) 

u XX ■ 



V. SHORT-RANGE SCATTERERS IN BORN 
APPROXIMATION 

The primary aim of this paper is to estimate the im- 

(2) 

portance of the contribution of a xx to the conductivity 

@xx • 

We have performed numerical calculation of these 
quantities for the realistic parameters of a bilayer system. 
The notably difficult problem is to model the electron 
scattering on impurities, as we have only limited knowl- 
edge about the nature of scatterers and their distribu- 
tion in the sample. For simplicity, we assume the short- 
range scattering of electrons on impurities randomly dis- 
tributed in both the left and right wells, as mentioned 
above. The concentration of impurities is assumed very 
low and the weak scattering on an individual impurity 
is treated in the Born approximation. The strength of 
the scattering on an individual impurity and the impu- 
rity concentrations are taken as adjustable parameters. 
Only the non-self-consistent version of the Born approxi- 
mation is used which leads to inaccurate description close 
to singularities in the field dependence of <r. 

In the course of calculations, we apply the standard 
procedure of averaging the resolvents and their products 
over all possible configurations of impurities. We assume 
(Vimp) = as usual. The averaged resolvent (G ± (i?)) 

is approximated by G (E), which satisfies the coupled 
equations 



G ± (E) 



1 



E - H - S ± ' 



(VimpG (E)Vi m p) 



(19) 



(20) 



where is the self-energy operator in the Born ap- 
proximation, i.e. determined to the second order in the 
perturbation Vi mp . 

The products (G ± v s G ± ), s = x, y, appear in the ex- 
pressions for conductivity. We approximate this expres- 
sions by the vertex functions K^ 1± , which must be found 



by solving the Bethe-Salpeter equations. In the Born 
approximation, the equations take the form 



(21) 



= G [V S + (Vi^K^Vimp) G 



where (VimpK^Vimp) are the vertex corrections to the 
velocity component v s . 

To proceed further, we accepted additional simplifica- 
tions motivated by the very low scattering rate in high- 
mobility samples. Namely, the real part of the self-energy 
S ± (a correction to the eigenenergies) is neglected and 
only the terms of the lowest order of the scattering rate 
are kept in all expressions. 

In this approximation, the imaginary part of the self- 
energy S 1 * 1 , denoted by T, is diagonal in the spectral rep- 
resentation of the Hamiltonian Hq and proportional to 
the densities of states, g a = m* /(-kTi 2 ), of the uncoupled 
electron layers, 



(22) 



7T7 L (1 ± A)g L + irj R (l =F A)g R . 



The parameters 7 Q represent the intra-well scattering 
rate and are given by the square of matrix elements of 
the single-impurity potential, multiplied by the concen- 
trations of impurities. Thus, Tf, a is determined as a 
weighted sum of scattering inside the left and right sub- 
systems with the weights 1 ± A. 

It is also worth to know how the resulting quantities de- 
pend on the scattering rate. While crjs and would di- 
verge in samples without impurities and the leading term 
of their power expansion is of the order 1/r, ais oc 1/r 
and cx 1/r, the expressions obtained as the response 
to Uh will stay finite in this limit, oc 1. Note that 
in spite of it a x 2 J oc 1/r as a xx is related to e' 1 ' through 
equation (Q) . 

The vertex corrections vanish for the velocity compo- 
nent v y , i.e. the product (G ± v y G ± ) can be approximated 

by G v y G . The matrix elements of v y are symmetric 
in k space, (a\v y (—k)\a) = —(a\v y (k)\a), their symme- 
try and the isotropy of the scattering by randomly dis- 
tributed short-range scatterers are responsible for van- 
ishing the vertex corrections in this case. 

This is not true for (G ± u 2; G' ± ). Due to the breaking of 
the time reversal symmetry by B _L x, (a\v x (—k)\a) ^ 
— {ot\v x (k)\a) and the vertex corrections have to be taken 
into account. We must solve two coupled Bethe-Salpeter 
equations for the vertex functions. Unfortunately, for our 
very simple model they are linearly dependent. Similar 
complication was found also in H, where the same model 
was used to describe the electron transport by the Boltz- 
mann equation. The second independent equation was 
supplied by the charge conservation rule SQr = —SQl- 
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VI. RESULTS AND DISCUSSION 

The parameters corresponding to the sample B studied 
experimentally by Simmons et al. El are used in our model 
calculations. The distance between wells d equals 135 A 
and the hopping integral t is supposed to be 0.9 meV. 
The density of electrons N e = 2.4 x 10 11 cm~ 2 yields the 
Fermi energy 4.29 meV. 

The applied in-plane magnetic field, qualitatively 
changes the topology of Fermi contours □ At zero field, 
the Fermi contours are two concentric circles. The larger 
circle corresponds to bonding states, the smaller one to 
antibonding states. The applied in-plane magnetic field 
shifts the centers of circles in the opposite directions in 
k space and gradually changes their shapes. The Fermi 
contour of the bonding subband acquires the shape of 
a "peanut" , the Fermi line of the antibonding subband 
evolves into a "lens" shape. The area of the lens is re- 
duced by the increasing in-plane field. At the critical 
field B c .i = 7.4 T the lens vanishes and the transition 
from two-component to one-component system occurs. 
The peanut splits at the higher critical field £? c .2 = 9.3 T 
into two parts, the left and right electron layers be- 
come decoupled and two disconnected ovals represent the 
Fermi contours of individual layers above that field. The 
changes in the connectivity of Fermi lines are reflected 
in the field dependence of the density of electronic states 
(DOS) as van Hove singularities at the critical fields. 
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FIG. 2. Calculated magnetic field dependence of the equi- 
librium diamagnetic current of electrons on the Fermi contour 
flowing in the right well. Except of the total current also the 
contributions from the bonding and antibonding subbands are 
shown: 6j x = Sjx.b + 5j x>a - 

As we have mentioned in the introduction, the in- 
complete balance of chiral diamagnetic currents flowing 
in the right and left wells is behind the origin of the 
electron transport in the x direction. In the linear re- 
sponse theory we implicitly assume that applied external 



fields are infinitesimally weak. Therefore, also the po- 
larization is weak and only the electrons with energies 
in a narrow interval around the Fermi energy contribute 
to the difference between the currents in the left and 
right wells. For this reason, Fig. || presents the calcu- 
lated equilibrium diamagnetic current 6jft. x flowing in 
the right well only for electrons on the Fermi contour. 
(The current of the opposite direction flows in the left 
well and Sjl.x = —&3r,x-) Together with the total cur- 
rent Sj x (= Sjn^x) also the contributions 5j x ,t and Sj x _ a 
from the bonding and antibonding subbands are shown. 
In our simplified model Sj x ^ and Sj Xja exactly cancels out 
for magnetic fields less than S c ,i- Above B c \ the current 
Sj x = 6j x .b decreases and changes sign before the mag- 
netic field reaches £? c .2- At that field Sj x diverges and for 
B > -B c . 2 returns slowly back to zero value expected for 
the decoupled layers. The sharp minima and maxima of 
the total Sj x reflect qualitative changes of the topology of 
Fermi contours, caused by the applied in-plane magnetic 
field, similarly as a field dependence of the DOS. 

The results for the conductivity tensor components a xx 
and a yy , calculated assuming zero temperature, are pre- 
sented in Fig. |^. The overall similarity of curves in Fig. || 
to 5j x is striking. The van Hove singularities at the crit- 
ical fields B c i and B c ^ seen on Sj x appear also on a xx 
and <jyy. Moreover, both components of the conductivity 
are small below B c \ where 8j x vanishes, are decreasing 
functions for B between B c \ and -B c ,2, and increase above 
-B c . 2 similarly as Sj x . But there are, of course, many dif- 
ferences which stem from the more detailed description of 
the scattering process and response to the external fields, 
as presented in previous paragraphs. 




B (T) 

FIG. 3. Magnetic field dependences of the conductivity 



components a xx and a y 



<7 (1) 
u xx 



is the conductivity calculated 



without the Hall correction (a x 2 J neglected). 



The dependence of the conductivity components a. 



(i) 
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and Oyy on the magnetic field coincides with that already 
obtained by solving the Boltzmann equation □. In this 
case the origin of the van Hove singularities is related 
rather to DOS than to Sj x . The difference between a2x 
and <7 yy originates in the Fermi contours anisotropy and 
the related anisotropy of the effective mass and the re- 

(2) 

laxation time. The novel "Hall" correction a xx to the 
conductivity a xx does not change the field dependence 
qualitatively, but strengthens the difference between a xx 
and a y y substantially. 

The question arises whether one can distinguish be- 
tween and (Txx experimentally. The possibility can 
be to study the inter-layer charge redistribution by mea- 
suring the change of the capacity between the bilayer 
system and the gate electrode on the top of the sample. 
It should be proportional to the current flowing through 
the sample perpendicularly to the magnetic field. The 
excess electric charge cumulated in the right layer, in- 
duced by the applied electric field in the x-direction, is 
shown in Fig. ^. The maximum effect is expected when 
only a bonding subband is occupied and the Fermi line 
has a single-connected-peanut shape. 
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APPENDIX: 

A Kubo formula evaluates the linear response of the 
density matrix f(t) to a perturbation of the form H p (t) = 
—AF(t), where A is a time-independent operator and 
F(t) is a function of time. The non-equilibrium ma- 
trix deviates from the Fermi-Dirac distribution function 
fo(H) by Sf(t). The expectation value of an operator B, 

SB(t) = Tr [BSf(t) ) , is related to the perturbation by 



SB(t) 



1 

ih 

x Tr 



dt'e et ' 



(fo(H) 



(Al) 
F(t + t'). 



Taking into account Hp , given by (fit), we get 




FIG. 4. Charging of the right layer when external electric 
field is oriented in the x direction. Note that Q is negative. 

To conclude: we have shown that the charge redistri- 
bution within a double-layer system gives rise to the Hall- 
like contribution to the conductivity which is so strong 
that it cannot be omitted. The experimental verification 
of this Hall effect is however nontrivial. Any potential 
contacts destroy the double-layer structure heavily and 
a direct measurement of the potential difference between 
layers can hardly be realized. We believe that the solu- 
tion is a non-invasive measurement of the charge redis- 
tribution from the capacity changes of gated structures. 



SB (t) = 









/ dt'e et ' 


(e iut> - l) 


+ 1 


I — CO 







(A2) 



xTr (f (H) 



e% Ht 'v x e-* Ht ' ,B 



where 1 was subtracted and added in square brackets of 
the first line of the expression. 

The integration by parts of the second term, corre- 
sponding to the added 1 in an above sum, converts ( [A2] ) 
to a form 



SB(t) = ^e Mt [° dt'e et ' ( e 1 ^' - 1 



(A3) 



+ 



x Tr [f Q {H) 



e% Ht 'v x e-% m ',B 



x,B 



To derive (A3), we employed v x = [x,H]/ih and the 
identity 



dt'ee^' Tr (/„(£) 



e* m 'xe-K Ht \B 



0, (A4) 



which can be proved replacing functions of Hamiltonian 
standing inside the trace according to prescription 



Z(H) 



dEZ(E)5(E — H), 



(A5) 



where Z denotes such a function, and performing the 
integration over the time variable. The result 



dE / dE'- 



ihe 



E' -E + ihe 
x fo(E) - f (E')] Tr S(E - H)x5{E' - H)B 



(A6) 



G 



is an expression which goes to zero when e — > 0. 

A similar approach is applied to the first part of 
(|A^), corresponding to the subtracted 1. After the 
time integration we introduce resolvents by G^(E) = 



( 



E — H ± iTiej , and, letting ui — > to get a static re- 
sponse, we obtain the Kubo formula in the form 



5B(t) = iheE x J dEfo(E) 



(A7) 




^dG+[E)_ 
dE 



dE 



+ l im ^ e -*Tr(7o(£) [Ml) 



The expression ([n]) for JQq, is derived by introducing 
Q a , defined by equation (||), into (A7) instead of £?, and 
taking into account that [x, Q a ] = 0. The conductivity 
formula (^) is obtained when B is replaced b y th e cur- 
rent operator (|l6|). In this case the last term of ( A7) gives 
a non-zero contribution, which exactly compensates the 
equilibrium gauge current. 

The expressions ( |l2| ) and ( p^| ) are derived in an anal- 
ogous way, taking into account the perturbation Hamil- 
tonian instead of Hp 1 * 1 . In fact, the derivation is 

simpler in this case. 
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